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A class of three dimensional classical lattice systems with macroscopic ground state degeneracies, 
most famously the spin ice system, are known to exhibit "Coulomb" phases wherein long wavelength 
correlations within the ground state manifold are described by an emergent Maxwell electrodynam- 
ics. We discuss a new example of this phenomenon — the four state Potts model on the pyrochlore 
lattice — where the long wavelength description now involves three independent gauge fields as we 
confirm via simulation. The excitations above the ground state manifold are bions, defects that 
are simultaneously charged under two of the three gauge fields, and exhibit an entropic interaction 
dictated by these charges. We also show that the distribution of flux loops shows a scaling with 
loop length and system size previously identified as characteristic of Coulomb phases. 



I. INTRODUCTION 

Much recent activity in condensed matter physics has 
involved the study of "topologically ordered" phases 
which characteristically exhibit emergent gauge fields 
and deconfined fractionalized excitations at low ener- 
gies. Canonical quantum examples of these are the var- 
ious fractional quantum Hall phasesii^ and much of the 
physics is present in the elegant classical physics of the 
spin ice compounds^. 

Gauge fields are intimately connected to local con- 
straints, as in the textbook example of Maxwell electro- 
magnetism wherein Gauss's law reflects a constraint at 
each point in space that must be obeyed by the dynamics. 
In a condensed matter setting, the analogous constraints 
arise as a low-energy effective truncation of the space of 
configurations; examples range from the dimcr configu- 
rations of short range RVB theory^ to the string nets of 
Levin and Wen^. 

A subclass of these constraints literally take the form 
of lattice versions of the familiar Gauss's law for abelian 
gauge fields, albeit with restricted microscopically real- 
izable values for the lattice electromagnetic fields. The 
introduction of an appropriate classical statistical me- 
chanics that consists of averaging over all allowed con- 
figurations with uniform weight leads to the so-called 
Coulomb phased with dipolar correlations, whose coarse- 
grained description realizes a Euclidean Maxwell theory. 
More elaborately, the introduction of a quantum dynam- 
ics in the constrained manifold can lead to a version of 
Maxwell electrodynamics coupled to electric charges and 
magnetic monopoles. An elegant point of intersection 
between the classical and quantum Coulomb phases is 
an appropriate Rokhsar-Kivelson point where the ground 
state wavcfunction is itself an equal amplitude superpo- 
sition of allowed configurations^i^. 

In this paper we expand the catalog of Coulomb 
phases. We study the antiferromagnetic four-state Potts 
model on the highly frustrated pyrochlore lattice and 
show that its ground state manifold exhibits correlations 



characterized by three abelian gauge fields. We find that 
the fundamental excitations/defects above this ground 
state manifold are charged under these gauge fields in an 
unusual way — they carry nonzero charges for two of the 
three gauge fields whence we refer to them as bions. (The 
term dyon is already reserved for particles charged un- 
der dual electric and magnetic fields whereas here both 
fields are of the same species.) In the classical setting, 
which is our primary interest in this paper, the import 
of the charge assignments is that it predicts the entropic 
force between different bions and more generally the free 
energy /entropy for any configuration of multiple bions. 
Our evidence for these assertions comes from a Monte 
Carlo simulation that agrees with the correlations pre- 
dicted by the triple Maxwell theory, and which yields 
statistics of flux-loops in the ground state manifold that 
have been previously suggested to be a sharp diagnostic 
of the Coulomb phased. 

Readers familiar with the existing literature on 
Coulomb phases will note that it is already know n^°d^ 
that classical 0{N) spins with a nearest neighbor anti- 
ferromagnetic interaction on the pyrochlore lattice 

(a) exhibit a Coulomb phase with one gauge field for N = 
1 (Ising) spins which is the case relevant to spin ice 
and indeed observed in experiments, 

(b) exhibit order by disorder for N ^ 2 and 

(c) exhibit a Coulomb phase with N gauge fields for 

> 3. 

It is interesting to situate the current work in this con- 
text. To this end imagine starting with 0(3) symmetric 
Heisenberg spins that live on the sphere (Fig. [5J:) whose 
ground state correlations arc governed by three indepen- 
dently fluctuating gauge fields. Excitations above this 
manifold are gapless and involve arbitrarily small charges 
under the gauge fields. If we restrict their range by gen- 
erating an easy axis (Fig. [2^) we return to the Ising case 
where the number of gauge fields is now down to one and 
the excitations are gapped and quantized. The import of 
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this current paper is that if we restrict the range to four 
symmetric points on the sphere (Fig. Oj) the number of 
gauge fields is unchanged although the excitations again 
become gapped and quantized. We believe that this re- 
duction can be extended to higher dimensions by con- 
sidering generalizations of kagome/pyrochlor involv- 
ing d + 1-simplices in d dimensions and starting with 
0(d) spins and restricting them to c?+ 1 state Potts con- 
figurations. 

We would be remiss if we did not note that this paper 
generalizes the early results of Kondev and Henlejsi^ from 
the two dimensional lattice known variously as the square 
lattice with crossings or planar pyrochlore, to three di- 
mensions. Readers who peruse the early paper will spot 
the family resemblance immediately. 

In the balance of the paper, we will set up the Potts 
model and its mapping to vector spins (Section II), map 
these in turn to a coarse grained description in terms of 
pseudo- magnetic flux/gauge fields and confirm the result- 
ing predictions for the correlations (Section III), discuss 
the bionic excitations (Section IV), study the statistics of 
loops (Section V) and conclude with some brief remarks 
(Section VI). 

II. THE MODEL 

The pyrochlore is a lattice of corner-sharing tetrahe- 
dra which can be constructed from the diamond lattice 
by placing a site at the midpoint of each bond (Fig. [1]). 
The result is a quadripartite structure, which may alter- 
natively be described as an fee lattice with a four-site 
basis. From the former construction, it is evident that 
the centers of the tetrahedra lie on the sites of the di- 
amond lattice: in other words, the dual lattice of the 
pyrochlore is the diamond lattice - a fact which we will 
make use of extensively below. We now turn to the Potts 
model which we will introduce from the perspective of 
the Heisenberg model as this will yield a vector spin rep- 
resentation of the Potts spins which will be central to 
this paper. 

The pyrochlore lattice is highly frustrated from the 
perspective of classical coUinear antifcrromagnctism: the 
nearest-neighbor classical Heisenberg antiferromagnet on 
this lattice has an extensive ground state degeneracy and 
remains a quantum paramagnet at all temperaturea^^. 
Since they will be relevant to us, we briefly summarize 
some details of the Coulomb phase for 0(3) (Heisenberg) 
spins on the pyrochlore. The canonical nearest-neighbor 
Heisenberg Hamiltonian 

i/-J^S,.S, (1) 

can be re-written, up to an overall constant, as 



where the sum in parenthesis runs over the four spins at 
the corners of each tetrahedron, and the outer sum runs 
over all tetrahedra in the lattice. Thus, the ground states 
are defined by spins satisfying local constraints: 

for each tetrahedron and each spin component a. 

As a result of these local constraints, each ground 
state can be mapped to a configuration of divergence-free 
fluxes, one for each spin component, on the dual diamond 
lattice. Upon coarse-graining, the entropic cost of fluctu- 
ations within the ground state manifold leads to an emer- 
gent 'electrodynamics' — with the coarse-grained fluxes 
playing the role of divergence-free lattice electromagnetic 
fields. The process yields asymptotically dipolar forms 
for spin (and field) correlation functions, a hallmark of 
the celebrated "Coulomb Phase" . 

We now consider applying a symmetry breaking poten- 
tial that restricts the Heisenberg spins to four symmet- 
rically situated points in spin space (Fig. [^b). The spin 
on each pyrochlore site must now belong to the following 
set of four spins pointing towards the corners of a regular 
tetrahedron in spin-space: 

Sa = (-1,1,1); Sb = (1,-1,1); 
Sc = (-1,-1,-1); Sz3 = (1,1,-1). (3) 

Observe that any two (different) spins in ([3]) make an 
angle of cos~^(— |) with one another, so that the dot- 
product of any two spins in the set is Sq • = 4 6ap — 1 
where a, /3 = A,B,C, or D. Thus the nearest neigh- 
bor interaction energy has the character of an antiferro- 
magnctic Potts interaction between four states: it prefers 
neighbors to be different but is indifferent to how that is 
achieved. Formally, the Hamiltonian (II} with the spins 
restricted to the set ([3]) is equivalent to the Hamiltonian 

Hp = jJ2S<r.,-r (4) 

where the Potts spins can be in one of four states: 
A,B,C or D. In essence we have mapped from Potts 
variables to a set of vector spins. The ground state con- 
straint ^ restricted to the set ^ is equivalent to the 
condition that each tetrahedron to contain all four Potts 
states. 

The Potts model has a discrete macroscopic ground 
state degeneracy — a remnant of the continuous macro- 
scopic degeneracy of the 0(3) modclM. To show this, a 
strict lower bound on the entropy can be obtained by us- 
ing the degeneracy of the three-state Potts model on the 
kagomoi^ii^ and by viewing the pyrochlore as alternat- 
ing layers of kagome and triangular planes. The result 
is the boundii fl > i (1.208 72)^/^^ corresponding to an 
entropy S/ksN > (1/2) log(1.208 72). A more direct es- 
timate is the Pauling entropyi^ for this system. For a 
given tetrahedron, 4! of the possible 4^ states are ground 
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states. Treating the tetrahedral constraints as indepen- 
dent gives a ground state degeneracy 



= 4 



N 



N/2 



where N is the number of spins and Ntet — N/2 is the 
number of tetrahedra. This corresponds to an entropy 
per spin S/ksN — (1/2) log(3/2) which is, interestingly, 
the same as the PauHng estimate for the entropy of spin 
ice^. As advertised in the introduction, the reduction 
from Heisenberg spins to Potts spins suggests that the 
latter system will still exhibit a Coulomb phase. We now 
turn to a precise formulation of this conjecture. 




III. FLUX FIELDS AND CORRELATIONS 

A. Flux Fields 

Our development of flux fields closely parallels the con- 
struction in the case of the 0{N) antiferromagne t ^ ^ . 
The essential idea is to map the spin variables to a flux 
field defined on the sites of the dual diamond lattice. The 
local ground state constraint - which still applies to 
the Potts spins as defined in ([3]) - is then mapped into 
a requirement that the flux configuration be divergence- 
free. 

We begin by defining bond vectors pointing from 
the even to the odd sublatticcs of the bipartite diamond 
lattice (i.e., from the centers of 'up' to 'down' tetrahedra 
on the pyrochlore), which take the form 
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(5) 



where the fee lattice constant has been chosen as a = 1. 
The spins live on the midpoints of the bonds; Fig.[l]illus- 
trates the geometry of the lattice and the bond vectors. 

Next, we define three vector flux fields on each bond, 
one for each spin component of the Potts spins repre- 
sented in ([3]): 



(6) 



where denotes the spin on bond u^, and a = 1,2,3 
labels the spin components. The flux field on a site of 
the diamond lattice is defined as the sum of the fields on 
the four tetrahedral bonds emerging from that site. 



B«(R) = V5^(RK, 



=1 



(7) 



where R is a diamond lattice vector, and k sums over 
the four tetrahedral sites (on pyrochlore) surrounding the 
diamond site (tetrahedron center). The mapping from 
spin to flux variables is invertible (see Appendix). 



FIG. 1. Geometry of tlie pyrochlore lattice. The centers of the 
up and down tetrahedra define the two sublattices of the dia- 
mond lattice, denoted by black and red dots. The pyrochlore 
lattice sites lie at the vertices of the tetrahedra. The vectors 
U-K define the bond vectors of diamond. 



From our definition of the Potts spins, we see that 
in any ground state, for each spin component S"", we have 
two "incoming" (-1-1) and two "outgoing" (—1) spins on 
each tetrahedron, i.e. each tetrahedron obeys a two-in, 
two-out 'ice rule' for each spin component. It follows 
from this, and our definition of the flux fields, that the 
local Potts constraint maps to a zero-divergence condi- 
tion for each of the B" fields, V • B" = 0. In an electro- 
dynamic representation it is appropriate to refer to these 
flux fields as 'magnetic' fields and then their sources will 
be monopoles — this is the nomenclature that is natural 
in the context of spin ice and is the one we will use heroic. 

Naively, the problem just looks like three copies of 
spin ice, one for each spin component. Of course, the 
three components are not independent and therefore we 
might expect some correlation between the three mag- 
netic fields. However, we will see shortly that our naive 
expectation is justified: at long distances, these cross- 
correlations vanish and the physics is described by three 
independent divergence-free 'Maxwell' fields. 



B. Coarse-grained Free Energy and Correlations 

Thus far, we have given a characterization of the 
ground state manifold in terms of divergence-free con- 
figurations of three magnetic fields. However, in order to 
compute correlations in the limit T — > 0, we need a more 
workable description of the ground state manifold which 
we will now obtain by coarse graining. 

Let us consider one of the three magnetic fields, say 
B^. If we were to fiip the direction of flux on one of the 
"in" bonds at a diamond site (say by switching spins Sa 
and Sb), the zero divergence condition would require us 
to also flip the direction of an "out" bond at the site. 
We can continue flipping spins in such fashion until we 
get either a closed loop of Sa, Sb spins, or a string of 
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Sa, Sb spins that extends across the entire system (for 
finite systems, the latter eventually closes through pe- 
riodic boundary conditions). Flipping spins that lie on 
closed loops leaves the net magnetic flux through the sys- 
tem unchanged, while flipping spins on spanning strings 
changes the net flux threading the system. Since the av- 
erage flux contributed by a closed loop is zero, systems 
with large numbers of closed "flippablc" loops will have 
a small net B^. On the other hand, a large and satu- 
rated net requires the field on almost every site of the 
dual lattice to point in the same direction and therefore 
the number of flippable loops is small - an intuitive pic- 
ture is that a the only flippable loops in a saturated field 
configuration are those that span the system. 

Thus far, we have only talked of lattice magnetic fields 
that live on the bonds of the pyrochlorc. To derive long- 
wavelength properties, we need to define smoothly vary- 
ing continuum fields. We do this by coarse-graining - the 
field B"(r) is defined as the average of the lattice fields 
B" in some neighborhood of r that is much larger than 
the lattice spacing but much smaller than the system 
size. The discrete constraint naturally translates into a 
divergence-free constraint for the coarse-grained fields. 

Microscopically, there are many more configurations 
consistent with a small net coarse-grained B^ rather than 
a large saturated B^; the same arguments obviously ap- 
ply to all three magnetic fields. Therefore configurations 
with small average fields are entropically favored. To low- 
est order, the (entirely cntropic) free energy as a function 
of the coarse-grained fields and consistent with symme- 
tries can be written as 



Ftot{B"{r)) 
_ 1 kT 



-TS 



Bi(r) 4- B2(r) 



B^ 



(8) 



where we've inserted a factor of the lattice-spacing a to 
make the stiffness, k, dimensionless. Since we are re- 
stricting our attention for the moment to ground state 
configurations, the coarse-grained fields still satisfy the 
zero-divergence constraint, V • B" = 0. The free energy 
([8]) coupled with the divergence-free constraint yields 
three copies of Maxwell electrodynamics in a standard 
fashion. Introducing three vector potentials A" to im- 
plement the constraints, we can rewrite the free energy 
as 

F=i^/dV^|VxA"(r)|^ 

^ •' Q = l 

We calculate the long distance correlators of the mag- 
netic fields, and the result is the dipolar form typical for 
Coulomb phases 

G^^(r)^(Br(r)Bf(0)) 
_ a ^ Srjrj - 

4:TTK 



(9) 



(a) 




where i,j refer to the x,y,z components of each 
magnetic-field. 



FIG. 2. Spin space representation of Ising, Potts and Heisen- 
berg spins (left to right). The red areas represent the ac- 
cessible portions of spin space - Ising spins are confined to 
lie on only two points, 4-state Potts spins occupy only four 
points, and Heisenberg spins can lie anywhere on the sphere. 
While Ising spins only require a single gauge field for their 
description, Heisenberg spins require three such fields. Some- 
what surprisingly Potts spins - which a priori would appear to 
have a much reduced symmetry compared to the Heisenberg 
case - also require three gauge fields. 



Finally, we reiterate that Ising spins, which occupy 
only two coUincar points on a sphere in spin-space, re- 
quire a single magnetic field B to describe their Coulomb 
phase. Classical 0(3) spins can lie anywhere on a sphere, 
and they require three fields B^, B^ and B'^ for their de- 
scription - one for each spin component. It is interesting 
that even though Potts spins are locked to just 4 points 
in spin space, they still require three magnetic fields for 
their complete description: the theory thus renormalizes 
at low temperatures to an effective 0(3) symmetry. Fig. [2] 
illustrates this idea. 



C. Monte Carlo Simulations 

In order to test our conjectured form ([9]) for the real- 
space correlation functions for the magnetic field, we turn 
now to a Monte Carlo study of the ground state manifold. 
We perform simulations on systems with L x L x L unit 
cells with L — 8, 16, 32 and periodic boundary conditions; 
since there are 4 sites per pyrochlore unit cell, this cor- 
responds to AL^ spins. The simulations were performed 
using a standard "worm-update" algorithm. In each up- 
date step we first identify a closed "worm" of alternating 
spin flavors (for instance AB . . . AB) , and then flip all the 
spins along the worm (i.e. interchange the spins A B). 
This move respects the Potts constraint since each bond 
in the worm is only part of a single tetrahedron, and 
exchanging spin-flavors on a bond still leaves a tetrahe- 
dron with all four Potts flavors. There are (2) = 6 types 
of worms, and the starting site for a worm and its type 
were chosen randomly for each update. It is instructive 
to think of this in the language of fluxes: interchanging 
A and B sites corresponds, via Eq. ([3]), to identifying a 
closed loop of type 1,2 and 3 fluxes and then reversing 
the first two of these. Reversing closed loops of fluxes 
clearly leaves the solenoidal constraints intact. 

We simulate M independent Markov chains, each with 
configurations along the chain {A'l and N were typi- 
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cally 100 and 10,000 respectively). To generate these, we 
begin with M independent "seed" ground state configu- 
rations and use the worm update described above to gen- 
erate the states along the chain. To ensure that succes- 
sive states along the chain were roughly independent, we 
perform several 30) such worm updates before record- 
ing a new configuration which is then added to the chain. 

In analyzing the data, we first obtain the average value 
of the correlation functions for each Markov chain, and 
then average these across all M chains. The error is esti- 
mated as the standard error of the single-chain averaged 
correlation function across the M independent chains. 

We compute the average correlation function G^''(r) = 

(Bf (x + r)_B^(x)) in two independent directions r for all 
36 combinations of a,(3,i,j. The vectors r were chosen 
as ri = ni (O, |, ^) and r2 = n2 (0, 0, 1) with ni,n2 & Z. 
These correspond to the fee lattice vectors ri = riiai and 
r2 = n2 (ai + a2 — aa ) where the a^ are fee basis vectors 
with lattice constant a = 1. 

The correlations fall off as 1/r^ consistent with the 
dipolar form ([9|). Figs. |3] and |4] show the representative 
correlators Gi^^(r), Gl.y{r), and G;J,^(r) multiplied by nf, 
in the two directions ri and r2 . The agreement with the 
dipolar form is best in the regime a r <^ L. The cross- 
correlators G^'^(r) for a ^ [3 vanish in all directions for 
all I, J, confirming the diagonal form for the effective free 
energy ([5]). We also checked that the ratios of correlations 
for different i, j asymptote to the values predicted by 
For example, G^?^/G?,^ —2 in the direction ri. 

Finally, we use correlation data to numerically esti- 
mate a value for the stiffness, n/a through ([9]). Cor- 
relations in the direction ri yielded an average stiffness 
Ki = 7.38 ± 0.46, while correlations in the direction Y2 
gave K,2 = 8. 13 ±0.93; the values in the two directions are 
equal within the margins of error. The lattice constant a 
is set to unity. 



IV. CHARGES, DEFECTS AND DIRAC 
STRINGS 
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FIG. 3. Monte Carlo data for three representative correlation 
functions G^^^iy^), G;5:^(ri) and Gyl(ri) in the direction ri = 
711 (0, |, ^) for lattice sizes L = 8, 16, 32. The correlators are 
multiplied by the cube of the distance nf , and the horizontal 
trends are consistent with the expected dipolar form (|9}. All 
cross-correlators C^^ for a ^ jS vanish (only Gj.^ displayed) 
confirming the diagonal form for the free energy ((8]). 
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While the effective free energies for the Heisenberg and 
Potts models are the same, the two differ in the nature 
of their excitations. Excitations above the ground state 
Potts manifold are gapped and quantized, with "bionic" 
defects that are charged under two of the three gauge 
fields. 

In the ground states, each site of the dual diamond 
lattice has two incoming and two outgoing fiuxes for each 
of the three magnetic fields. We can create defects by 
violating the zero-divergence constraint at a dilute set of 
points in the lattice. Such defects are "charged" under 
the different fields, with positive (negative) charges equal 
to the net outgoing (incoming) fluxes of each type at the 
defect. This is the usual charge in the sense of Gauss's 
law: each such charge represents a source of magnetic 
flux, and a violation of the divergence-free constraint for 
at least one of the fields. 



FIG. 4. Monte Carlo data for three representative correlation 
functions G]^^{r^), Gi^(r2) and Gyz(r2) in the direction r2 = 
712 (0, 0, 1) for lattice sizes L = 8, 16, 32. Once again, the 
correlators are multiplied by the cube of the distance 712 , and 
the trends are consistent with the expected dipolar form ([9}. 
Note that the correlations are weaker in the direction r2 as 
compared to ri because same values of n for the two cases 
correspond to larger physical distances r2. 



It is convenient to first catalog defects in the Potts lan- 
guage: a defect arises when the four spins surrounding 
a dual lattice site violate the Potts rule. The simplest 
defects are those in which one spin flavor is repeated on 
a tetrahedron; for instance, we can have a defect tetra- 
hedron with spins Sb, Sb, Sc and Sd- There are twelve 
such defect tetrahedrons: we have four choices for the 



6 



23 



Type of Defect 



Charges and Flux Lines 




FIG. 5. Charges of the twelve types of bions (black dots) 
shown in Qi,Q2,Q3 space, where Qi represents the charge 
under the field B*. Each bion is charged under two gauge 
fields. The twelve charges naturally map to six oriented edges 
of a tetrahedron. 



spin flavor that gets repeated {Sb in our example), and 
three choices for the flavor that the repeated spin replaces 
{Sa in the example). 

Each of these defects have different charges under the 
three gauge fields. Looking at the spins in ([3]), we see 
that an "up" tetrahedron with Sb, Sb, Sc, Sd has spin 
components S*^ = (1, 1, -1, 1), = (-1,-1,-1, 1) and 
= (1,1,— 1,-1). Thus, its charges under the three 
magnetic fields B^, and B'^ are Qi = +2, Q2 = —2 
and Qs = respectively. All twelve defects have a 
similar structure, in that they are doubly charged un- 
der two of the three magnetic fields - hence the name 
bions. Table |T] catalogues the charges of the different 
bions. (These charges arc reversed for corresponding de- 
fects on "down" tetrahedra, since the sense of "in" and 
"out" flux is reversed). Figure [5] depicts the 12 possible 
bions in Qi, (32, space. 

Charge conservation demands that the defects are al- 
ways created in oppositely charged pairs. One way to do 
this is to imagine creating a pair of bions by exchang- 
ing spins along a "Dirac string" containing two flavors 
of spin. For instance, we create a BBCD defect by re- 
placing an A spin with a B spin on an "up" tetrahedron. 
This creates a second, oppositely charged defect on the 
adjoining "down" tetrahedron. The second defect can be 
moved away from the first by continually flipping a string 
of A, B type spins. The second defect is of type AACD 
when it lies on another "up" tetrahedron; we may ver- 
ify that the tetrahedron AACD carries opposite charge 
Qi = —2, Q2 = -1-2 and = 0. We can think of the 
Dirac string as a flux tube carrying two flavors of flux 
(B^ and B^ in our example) that connects bions that are 
oppositely charged under two magnetic fields. Figure [6] 
shows a Dirac string connecting two bions. 

Additionally, we can also imagine creating composite 
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A A B D 



C C B D 
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: -2 
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: +2 
: -2 



: -2 


: -2 
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Connected by a Dirac string of 
> type B D 



Connected by a Dirac string of 
type B C 



Connected by a Dirac string of 
type C D 



TABLE I. Catalog of defects and charges. Each defect is 
charged under two gauge fields - hence they are called bions. 



defects by adding two or more of the twelve fundamental 
bions. The composite charges form an fee lattice which 
is a natural extension of Figure [5] There is an intuitive 
geometric picture for understanding the charge structure 
of the bions. The four spins point to the four corners of a 
tetrahedron in spin-space (and all spin-components sum 
to zero); replacing Sa with Sb to create a defect gives 
a vector of charges Q = — S^i -I- Ss under the different 
gauge fields. The charge Q thus corresponds to an edge 
of the spin-space tetrahedron. In this way, all twelve 
fundamental bions can be mapped to six oriented edges 
of a tetrahedron in spin-space. 

The bions act as sources of magnetic flux and expe- 
rience a Coulomb force under each of the three mag- 




FIG. 6. Dirac string connecting two bions. The colors orange , blue, green and red represent the spin flavors Sa, Sb, Sc and 
Sd respectively. The three arrows on each site represent the three magnetic fields (black), (purple), and B^ (gray), 
(a) Ground state conflguration of the 4-state Potts model with no defects. Each tetrahedron has all four spin flavors, and all 
three magnetic flelds obey two-in-two-out rules at each tetrahedron, (b) Switching Sa (orange) to Sb (blue) on the bottom 
tetrahedron creates (BBCD) type bions on two adjoining tetrahedra (red and blue spheres). The magnetic fields are no longer 
divergence-free, and the bion on the up tetrahedron has charges Qi = +2, Q2 — —2, Q3 = (opposite charges for the bion 
on the down tetrahedron), (c) The bions move apart by flipping a trail of Sa, Sb spins. The Dirac string acts as a flux tube 
carrying B^ and B^ type magnetic fluxes between the two bions. 



netic fields. The origin of this force is purely entropic 
in nature. We imagine a finite number iV^ of bions scat- 
tered throughout the lattice at positions ri, r2, • • • tat^. 
We can compute the partition function Z by integrating 
over all configurations of magnetic fields consistent with 
the distribution of bions. This gives the free energy of the 
sources in accordance with Z/Zq = e^^"^*^"^ , where Zq is 
the partition-function in the absence of bions. Explicitly 
this yields: 



F 

^ IT. 



y — 

^ Ann 

i=l 
j<i 



aT [QiQi + QlQi + QlQi 



E' 

i=l 



aT (Q^ 



inn 



(10) 



where i,j sum over all pairs of bions and represents 
the charge of the ith bion under the field B". The first 
term in (|10p represents the Coulomb interaction energy 
of pairs of bions separated by distance jr; — rj|, while 
the second term is the free energy of individual, isolated 
bions arising from self-interaction in the field theory. The 
self-interaction term has an ultraviolet ambiguity, rep- 
resented by the unknown constant a, which (naturally) 
cannot be fixed by the coarse-grained free energy alone. 
Also note that ([TU)) predicts that the interaction en- 



ergy between a pair of bions is sensitive to their type. 
For example, bions of type AABC and DDBC (equally 
and oppositely charged under Bi and B3) interact more 
strongly than say AABC and AABD which arc charged 
under different gauge fields. 

This is a good place to briefly comment on the finite- 
temperature properties of the Potts spins - in particular 
the form of the correlation function Eq. We know 
that the dipolar form of the correlation function is de- 
rived from the divergence-free constraint on the mag- 
netic fields. This constraint is exactly satisfied at T = 
and gradually weakened as the temperature T is in- 
creased. Heuristically, we might expect the correlation 
to be dipolar up to some (temperature dependent) cor- 
relation length ^(r), and decay exponentially on length 
scales longer than ^. It is easily seen^S that for Ising 
spins, the creation energy of gapped ice-rule violating 
defects (monopoles) yields ^ ^ e^"^^"^^ . On the other 
hand, for Heisenberg spins the gapless excitations yield 
much softer violations of the divergence-free rule and, 
correspondingly, a much shorter finite-temperature cor- 
relation length given by^i ^ ^ Xj \/T . For gapped Potts 
spins, the creation energy of a Bionic defect is 4 J which 
gives ^(T) ~ e^'^/'^'^; the gap to excitations helps preserve 
the dipolar form of the correlations to higher tempera- 
tures. As is typical, the gap also manifests itself in the 
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exponential low-temperature decay of various thermody- 
namic quantities but these are not the focus of this paper. 

V. WORM LENGTH DISTRIBUTIONS 

Coulomb phases come with a natural incipient loop 
structure — absent defects one can define closed lines of 
flux thanks to the underlying conservation laws. This 
makes the statistics of the loops worthy of interest. In- 
deed Jaubert, Haque and Moessncr (JHM)^ have stud- 
ied loop statistics for the ground state manifold of spin 
ice and found a characteristic scaling of their probabil- 
ity distribution which they have related to the proper- 
ties of random walks in three dimensions. This suggests 
that this scaling might be more generally associated with 
Coulomb phases and we will investigate and verify that 
possibility here. 

Specifically, JHM have studied the distribution of 
worm lengths for spin ice which is easily done by keeping 
track of the worms used to update configurations in the 
Monte Carlo. As in our problem, worms in spin ice are 
closed strings of alternating spin flavors but now with 
the feature that while there is only one species of worm, 
at each step there is a binary choice that must be made 
randomly^^. For these JHM found a characteristic scal- 
ing for probability distribution of worm lengths p{£), 

m = ia / (^) (11) 

where l refers to the worm length and L is the linear 
dimension of the system size. This scaling form unifies 
two populations of worms: short worms whose probabil- 
ity scales as p{t) ^ and long winding worms (that 
close after spanning the system through periodic bound- 
ary conditions) for which p{l) ~ . 

We have investigated analogous worm-length distribu- 
tions in the Potts model. It should be noted that the 
Potts model has six species of worms (each worm has 
only two Potts spin flavors and (2) = 6). However, by 
symmetry, all six worm types have identical distributions 
in the Coulomb phase. Fig. [7] shows the worm-length 
distributions obtained for system sizes L = 8, 16, and 32 
plotted in scaled variables. Leaving aside the deviations 
at very small and very large loops sizes we see that that 
the loop distribution indeed obeys the scaling form pT|) 
which is then independent evidence that the Potts model 
is in a Coulomb phase. We direct the reader to Ref. H for 
a rationalization of this scaling in terms of the properties 
of random walks. 



VI. CONCLUDING REMARKS 

To summarize, we have shown that the classical anti- 
ferromagnetic four-state Potts model on the pyrochlore 
lattice is in a Coulomb phase described by three emer- 
gent gauge fields as T —s- 0. It is instructive to view 




10'^ 10^ 10' 10° 10^ 10^ 



FIG. 7. Probability distribution of worm- lengths I for system- 
sizes 1/ = 8, 16, 32. We clearly see two regimes - an scal- 
ing for short loops, and an I independent scaling for winding 
loops. 



the Potts model as arising from a symmetry breaking 
potential that restricts 0(3) spins to just four points in 
spin space. Nevertheless, we have shown that the long- 
wavelength effective free energies of the 0(3) and Potts 
models are identical. 

An important point of difference between the Heisen- 
berg and Potts models lies in the nature of excitations 
above the ground state manifold. While the Heisenberg 
model involves gapless excitations, the Potts model has 
gapped excitations with a novel charge structure. We 
find twelve types of "bionic" defects, each charged under 
two of the three gauge fields. The charges are deconfined 
and can be connected by "worm-like" flux tubes of al- 
ternating spin flavor. We computed probability distribu- 
tions for lengths of closed worms, and found scaling laws 
in accordance with previous diagnostics of the Coulomb 
phase. 

The evident next step in this program is to incorporate 
quantum dynamics in a quantum Potts model exhibiting 
a Rokhsar-Kivelson point. We expect to report the de- 
tails of this construction and a study of the excitations 
and phase diagram in a future publications^. 
Note added. - Upon conclusion of this work, we became 
aware of overlapping results of a study of the same 
problem^i. 
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Appendix A: Spin-Spin Correlation Functions 



as 



It is useful to have a referenee for eonverting the flux- 
field correlation functions ([9|) to spin-spin correlation 
functions for the Potts spins living on the four sublat- 
tices of pyrochlore. We have three flux fields B^, B^, 
B'^ each with three components, for a total of nine flux 
components; these are labeled for a = 1,2,3 and 
I = x, y, z in One the other hand we have four Potts 
spins on every diamond site, each with three components, 
for a total of twelve spin components; we label these 
for a = 1, 2, 3 and k — 1,2, 3, 4. However, the constraint 
equations on the Potts spins ([2]) ensure that there are 
only nine independent spin components, thereby permit- 
ting an invertible mapping from flux-fields to spins. 

As an example, we use the definitions of the fiux fields 
on the diamond sites ([7]) and the definitions of the bond 
vectors ([5]), to explicitly write the x,y,z components 
of 

B\v) - 5f (r)ui + 5f (r)u2 + (r)u3 + (r)u4 



Bl 



^3 



Sa 



( QX QX QX I QX\ 



(5f 



nx QX 

O2 



where we have dropped the explicit dependence on r to 
simplify our notation. Finally we impose the spin con- 
straint equations ^ to eliminate 5" and we get 

Bl = \ {-SI - SI) 

Bl = I {-Si - SI) 

Bl = \ {Sf + Si) . 

A similar exercise can be carried out for the B^ and B^ 
spins to give the following matrix of transformations: 



Bl 
Bl 
Bl 
Bl 
Bl 
Bl 
Bl 

\bV 



V 



-i 



o\ 




(si\ 







Sf 







SI 







5f 







SI 







SI 


1 

2 




QX 


1 

2 




SI 

\sV 


0^ 





We can invert the relation above to obtain the spins in terms of the fiux fields as follows 

(B}\ 



(Sl\ 

s\ 

SI 
SI 



QX 

qV 

'-'3 

\sV 



/-I 




1 




-1 



V 





1 




1 




-1/ 



Bl 
Bl 
Bl 
Bl 
Bl 
Bl 

Bl 
\Bl) 



V 



These relations allow us to express spin-spin correlation 
functions as simple linear combinations of the -Bf corre- 
lations. 



Appendix B: Symmetries and Stability of the Free 
Energy 

We have claimed in the bulk of the paper that the 
long-wavelength coarse-grained free energy is insensitive 
to breaking Heisenberg (0(3)) symmetry down to a Potts 
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symmetry. In this section, wc formalize this claim using 
a renormalization-group argument. We will show that 
the lowest order terms that can be added to the free 
energy (and are allowed by symmetry considerations) are 
irrelevant in an RG sense. 

Consider a given ground state configuration of Potts 
spins on the pyrochlore lattice. Permuting the spins (say 
by exchanging spins of type Sa and 5*^ on every tetra- 
hedron) gives another ground state configuration. Each 
such "internal" symmetry of the microscopic spin de- 
grees of freedom should translate into a symmetry of the 
coarse-grained B fields^ in accordance with Appendix 
A. In fact, six permutation group elements are required 
to generate all 4! = 24 permutations of the Potts spins; 
these map to the following symmetries of the free energy: 

/(Bi,B2,B3)-/(B2,Bi,B3) 
= /(B3,B2,Bi) 
= /(Bi,B3,B2) 
= /(-Bi,-B2,B3) 
= /(-Bi,B2,-B3) 
= /(Bi,-B2,-B3) (Bl) 

An example best illustrates how one arrives at the sym- 
metries in Eq. (|Bip . A given tetrahedron has magnetic 
fields defined by 

= Via + Ufc - Uc + Ud 

B2 =Slna + si Ub + Sln, + Si xia 

= +VLa - Ub - Uc + Ud 

B3 ^ S\Ua+ S'g Ub + 55 Uc + 5|, Ud 
+Ua + Ub - Uc - Ud 

where {a,b,c,d} £ Permutationjl, 2, 3, 4} label the u 
sublattice vectors on which the spins {Sa, Sb, Sc, Sd} 
live. Then, exchanging spins Sa and Sb leads to the 
modified fields 

B; = 5| u„ + S% Ub + 5g Uc + Sf) Ud 

= +Ua - Ub - Uc + Ud 
= B2 

B'2 =SlUa+ S\ Ub + Slu, + Si Ud 
= Ua + Ub - Uc + Ud 
= Bi 

B^, =S'BUa+ S\ Ub + 5J U, + 5|, Ud 
= -f Ua + Ub - Uc - Ud 
= B3 



with Bl and B2 exchanged. The same analysis carries 
through for all tetrahedra, and exchanging Sa and Sb 
everywhere on the lattice is equivalent to exchanging the 
coarse grained fields Bi and B2. This gives the first 
of the symmetries listed in Eq. (|B1[) : the others can be 
derived in an analogous manner. 

Let us pause to consider the implications of Eq. (jBl[) . 
The first three equations require a symmetry under ex- 
changing any of the B^ fields; this justifies having the 
same coefficient in front of all three quadratic, diagonal 
terms in the free energy Eq. ([5|). The last three equa- 
tions forbid quadratic terms that are off-diagonal in the 
fields. For example, a term like B1B2 is not a symmetry 
under {61,62,63} {-61,62,-63}. At this point 
quadratic terms like B'fB\ can still be added, though we 
will soon show that these are forbidden by spatial sym- 
metries. 

Having considered transformations of the "internal" 
spin degrees of freedom, we now turn to lattice trans- 
formations. The space group FdZm of the pyrochlore 
lattice consists of the 24 element tetrahedral point group 
43m and a further 24 non-symmorphic elements, involv- 
ing a combination of rotations or reflections with trans- 
lation by half a lattice vector. For the pyrochlore dressed 
with Potts spins, the space group elements transform one 
ground state configuration into another. 

The 6q, , are lattice vector fields whose transformation 
under the space group elements (like rotations R) derives 
from the transformation of the lattice bond vectors u. 
For example, the field 61 transforms as: 

4 

61(r) = 5]5;^(r)u^, 

K=l 

4 

^5]5^(i?r)i?^^u^, 

Since the free energy involves an integral over all space, 
the change in the spatial location of the fields (r — > i? r) 
can be undone by a simple change of integration vari- 
ables; what matters is the transformation of the vector 
indices of the fields. Microscopically, the transformation 
of the vector indices derives entirely from a permutation 
of sublattice indices i.e. a spin belonging to sublattice 1 
of a tetrahedron at location r gets rotated to, say, sublat- 
tice 3 of the tetrahedron at location R r. In this way, all 
we're concerned about is the action of the space group 
elements in permuting sublattice indices. The 4! per- 
mutations of the sublattice indices lead to the additional 
symmetries: 

f{B:,By,B^j^fiBy,Bi,B^j 

= f{Bl,Bl,Bl) 

= f{Bl,Bl,Bl) 

= f{-Bl,-Bl,Bl) 

= f{-Bl^Bl^-Bl) 

= fiBZ,~Bl,~B^J (B2) 
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where a = 1,2,3 labels the type of magnetic field and 
we have used a compressed notation to label the nine 
arguments of the free energy function. 

We should emphasize that permuting sublattices is 
very different from permuting spins. In the latter case, 
we exchange spins of types Sa and S b regardless of the 
sublattices on which they lie; in the former, we exchange 
the spins living on sublattice 1 and 2 regardless of their 
type. As shown by Eqs. (jBip . (|B2[) . exchanging spins 
leads to symmetries under exchanging different types of 
B fields, while exchanging sublattices leads to symme- 
tries under exchanging different spatial components of a 
given type of field. 

As before, let us consider an example to understand the 
symmetries listed in Eq. (jB2p . Fix the center of one tetra- 
hedron as the origin and rotate the lattice by tt about the 
axis (0,0,z). (This is the Cz element of the tetrahedral 
point group). This axis bisects the edges Ui2 = (ui — U2) 
and U34 of the tetrahedron at the origin. The rotation by 
TT does a dual exchange of sublattice indices 1 2 and 
3 4 on each tetrahedron (in addition to rotating the 
tetrahedron's center). A given tetrahedron has Bi fields 
defined by 



^3 U3 



= Q.2b{-Sf + -31 + SI) 



Slul 



Bl 



SI ul + S^ ul + Slul + SI ul 
Q.2b{Sl~ S'^- Sl + SI) 



Bl = 



SI ul + 
0.25(5^ 



5"? 



S^ 1*3 



'^2 '-'3 



St) 



where {Si, S2, S3, S4} are the spins living on sublattices 
{ui, U2, U3, U4} respectively as in Appendix A. Now, ex- 
change Ul ^ U2 and U3 ^ U4. The transformed field 



equations are: 



{BfY ^Sfu^^+S^uf 



Si 



= 0.25{+S^ -S^ + SI- S^) 
= -B^ 



{Bfy 



Sful 
0.25(- 



f 5f 



SI 



SI u\ 



SI 



Sf) 



-Bl 



[Bl)' = + ^2" ul + S'iul + Sl'. 



0.25(5? 
Bl 



5? 



QX 

•^3 




The transformation has flipped the sign of the first two 
components of Bi corresponding to the fourth symmetry 
in Eq. (|B2p . Of course, exactly the same transformations 
carry through for B2 and B3. Also note that in 3D, the 
matrix for rotating by tt about the z axis looks like 



whose action is also to flip the first two vector indices of 
the fields it acts on. 

Finally, armed with the symmetries in Eqs. ()Bip . (|B2|) 
it is easy to see that the simplest terms we can add to 
the quadratic free energy defined in Eq. are cubic in 
the fields, and symmetric with respect to exchanging the 
types of fields and the x, y, z components of each field: 

F = Fquad + j (fr{BlBlBl + BlBlBl+^erumta.tions). 

(B3) 

The added terms are cubic in the gradient of the vector 
potential and thus irrelevant under a renormalization- 
group analysis for determining the long-wavelength cor- 
relations of the fields. This confirms the stability of the 
quadratic, diagonal free energy used in the bulk of our 
analysis. In future work, it would be interesting to explic- 
itly derive the perturbative corrections to the correlations 
stemming from Eq. (|B3[) . 
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